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^S) ' Abstract 

We investigate the dynamics of a Bianchi I brane Universe in the pres- 
ence of a nonlocal anisotropic stress Vij,i, proportional to a " dark energy" 
U. Using this ansatz for the case U > we prove that if a matter on a 
' brane satisfies the equation of state p ~ {y — l)p with 7 < 4/3 then all 

, such models isotropize. For 7 > 4/3 anisotropic future asymptotic states 

' are found. We also describe the past asymptotic regimes for this model. 

o 

O ■ 1 Introduction 

The idea that our space has more than three dimensions was attractive dur- 
ing many years, mostly in the scenario of compactification. Recently a new 
bX). paradigm appeared. In this new approach the matter is trapped in three- 

dimensional brane embedded in a multi-dimensional space (bulk) In the 
case of one extra dimension it is shown that isotropic cosmological evolution 
being unusual in early stages of the brane Universe ||] tends to the standard 
cosmological scenario in the late time 

The field equations on the brane are jj, |j 

Gf_ii, = K^Tfw + li'^S^i, — E^v (1) 

Here k is the fundamental 5-dimensional gravitational constant, n is the effec- 
tive 4-dimensional gravitation constant on the brane. Bulk corrections to the 
Einstein equations on the brane are of two forms: there are quadratic energy- 
momentum corrections via the tensor S^^, and nonlocal effects from the free 
gravitational field in the bulk S^^- The quadratic corrections are significant 
only in the early stages of cosmological evolution (see the next section). The 
nonlocal corrections can be decomposed into a scalar (nonlocal energy den- 
sity), a vector (nonlocal energy flux) and a tensor (nonlocal anisotropic stress) 
parts 1^. On an isotropic brane there is only a scalar part U which decays 
during the expansion of the Universe a,sU ^ a^* where a is the scale factor. 
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However, if we remove the suggestion of isotropy (for example, in order to 
study the problem of isotropization of the brane Universe), then the tensor part 
of nonlocal contribution P^j^ should be taken into account. The main obstacle to 
do this is the absence of an evolution equation for the anisotropic stress V^v 
^. Its correct form should be derived from a full 5-dimensional anisotropic 
metric which is still unknown. In this situation we have to choose some ansatz 
for Vfi^v In first papers on Bianchi I brane either all nonlocal corrections Q 
or its tensor part ||, § have been neglected. This approach gives some 
interesting results such as isotropic character of the initial singularity Q and 
the possible existence of an intermediate anisotropic stage However, whether 
these results remain valid when the contribution from V^v becomes important 
is still an open problem. 

Recently Barrow and Maartens have proposed an ansatz for the nonlocal 
stress which generalizes known anisotropic sources in General Relativity |pl| . 
They have studied the dynamics of isotropization in the limit of small anisotropy. 
In the present paper we use this ansatz for describing the global properties of 
the dynamics. We use the techniques of dynamical systems and the expansion 
normalized variables (for the description of this formalism see [|l^) to find the 
past and the future attractors of our system. The transient behavior is also cov- 
ered by studying the future asymptotic states for a pure nonstandard (i.e. with 
Sfj^y and without T^^, terms in (1)) brane dynamics. Of course, in this frame- 
work we can not find regimes originated from an interplay between standard 
(proportional to T^^) and nonstandard (proportional to S^i,) source terms. Up 
to now several interesting regimes of this type have been found (for example, an 
analog of the Einstein static Universe and a stable oscillationary regime Q ) . 

The paper is organized as follows: in Sec. 2. we describe the dimensionless 
variables and use them to write down the system of differential equations in 
a convenient form. In Sec. 3 we list stable points of this system neglecting the 
square energy-momentum correction and describe the future asymptotic states 
of the Bianchi I brane model. In Sec. 4 the stable points of the system without 
linear energy-momentum contribution are listed and possible transient behavior 
is obtained. We also discuss the past asymptotic states of our model. Sec. 5 
provides a brief summary of the results obtained. 

2 Expansion normalized variables and equation 
of motion 

We will use an orthonormal frame of vector fields Ca such that cq ~ dt- In the 
Bianchi I case we have the following commutation relations for the vector basis 
(Latin indices run from to 3 and Greek indices from 1 to 3): 

[ea, ef,] = C^jBc 

where 
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Here fl" is the angular velocity of the spatial triad with respect to a Fermi- 
propagated frame and Qfj,^, can be decomposed into a volume expansion rate G 
and a shear cr^^: 
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If we introduce the metric for the Bianchi I model in the standard form 



ds-" = -dt^ + af{t){dxy, (2) 

then the volume expansion rate in terms of the main scale factor a = (010203)^/^ 
is = 3H = 3^, where H is the main Hubble parameter. 

For these variables we can write the Raychaudhuri equation 



the Gauss-Codazzi equations 



CT^i. + Bcr^^ = -^"P^i., (4) 

Hj a 

--0 +0-'^ (Tp,.+2K p = -K — (5) 

and the conservation equations 

p + e(p+p) = o, (6) 
+ ^ew + (T^'^p^, = 0, (7) 

= 0. (8) 

Derivation of these equations from Eq.(l) and the explicit form of S*^^ see in 
Here A is the brane tension. Quadratic energy-momentum corrections (first 
terms in RHS of Eqs. (3) and (5)) become significant if a matter density p on 
the brane is greater than A. The current experimental limit on the brane tension 
is A > \m{GeVf 0. 

Due to presence of the nonlocal stress 'P^lv this system of equations is not 
closed. Therefore, it is necessary to specify the form of 7^^^ using some addi- 
tional suggestion. In the present paper we deal with the proposal of Barrow and 
Maartens [|ll]. They have suggested that the nonlocal stress is proportional to 
the nonlocal energy density 

= ID^^JA, (9) 
D^jy = diag{ci, C2, -(ci -I- C2)}, 
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where ci, C2 are some constant. 

It is worth to notice that choosing the appropriate coordinate system we can 
diagonaUze either Vfj,i, or cr^i,. We choose to make 7^^^ diagonal. 

Nevertheless, our system of equations is not completely determined yet. We 
must fix the angular velocity fl" of the spatial triad associated with V^p. After 
that we can write the derivatives of the shear tensor as 

In the GR description of cosmological anisotropic sources (for example, the 
homogeneous magnetic field |l5|] ) the angular velocity components are con- 
strained by the corresponding equations of motion of the anisotropic fluid (for 
example, the Maxwell equations in the magnetic field case). As we have no 
equation for T^^^y, we should choose an ansatz for Q". In the present paper we 
consider the simplest possible form fl" — (in the end of the next section one 
important statement not depending on fi" is proved). 

After introduction of the deceleration parameter 

and dimensionless variables 

E_ 

E23 
E13 

U 

we may rewrite our system in a simpler form. 

We have also auxiliary evolution equation for the matter density. There 
are two ways for description of the quadratic energy-momentum terms. It is 
possible either to incorporate them into a single evolution equation for the 
matter as it have been done in pfl or to consider separately as some kind 
of a particular effective matter 0, We choose the latter case and write two 
evolution equations - one for " standard" matter and the other for " nonstandard" 
matter. In order to do this we introduce dimensionless variables 



aa 



(10) 



-((T22 + o-33)/e, 

i\/3(cr22 - Cr33)/e, 

\/3(7i2/e, 

V3c723/e, 

V3ai3/e, 

mi 



(11) 



Note that these values are not independent, they connected to each other by 

n^/n^ = p/2X. 

Finally, we introduce a new time variable t such that ^ = ^- Now our 
system becomes (prime denotes a derivative with respect to r) 

/3E_(ci - C2)], 



(13) 



u' 


= 2C/[(<7-l)-3E+(cH 


-C2; 




= 3U{ci+C2) + Y.+ {q- 


2), 




= V3J7(ci-C2) + S-(? 


- 2 




= (9-2)Si2, 




^13 


- (9-2)213, 




y 

^23 


= (9 - 2)S23 





with the constraint 



the equation for q 



(14) 



g = 2E2 + i(37-2)f7^ + (37-l)f}A + t/, (15) 
and the evolution equations for matter density 

n'^ ^ [2{i + q) - (16) 

n'^ ^ [2(1 + q) - 6j]nx. (17) 

Here 

= + + E?2 + + E?3. (18) 

There is an important difference between GR and brane problems. In com- 
mon GR formulation of the problem the energy density of an anisotropic matter 
is positive, though the nonlocal energy density on the brane can be either posi- 
tive or negative and, in principal, can change sign in the cosmological evolution. 
However, ansatz (9), as we can see from the first equation of system (13) prevents 
U from the change of the sign. The case of negative energy of an anisotropic 
matter has not been investigated previously in the GR formulation due to its 
peculiarity from the physical point of view, but in the brane scenario it becomes 
important. The main feature which arises from U < is that the flat Bianchi 
I Universe can recollapse without reaching any future nonsingular asymptotic 
regime. From the mathematical point of view the negativity of U makes our 
phase space noncompact (the phase space, specially constracted to be compact 
in the U < case see in [|j). So, from this point we assume U > leaving the 
opposite case to a future work. Now our phase space becomes compact with 
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S+, E_, Si2, Si3, S23 varying in the range [—1, 1] and O^, f^A, U varying in the 
range [0, 1]. 

In the next section we start to analyze the system (13) - (17) neglecting 
the "nonstandard" matter. All the future attractors for this system are also 
the future attractors for the initial general problem, because for t ^ 00 we can 
neglect a nonstandard matter on the brane in comparison with a standard one. 

3 Equilibrium points and future evolution 

Now we list equilibrium points for the system (13) - (17) in the case Cl\ = 0. 
It is convenient to use the notations c+ = ci + C2 and c_ = ci — C2- We also 
have found the corresponding eigenvalues Aj of the linearization of system (13) 
about the equilibrium points and have used the condition Aj < to study the 
stability of the equilibrium points in the future direction. 

1. The isotropic Universe 

U = = S_ = S12 = S13 = S23 = 

g=^7-l f2^ = l 0<7<2 

Stable for 7 < 4/3 independently on c_|_,c_. 

2. Equilibrium point 

(2-7X37-4) _ c+(37-4) 
124 +4c^ + 64 + 2c^ 

_ V3c-(37-4) 37-4 
18cl + 6cl 184 + 6c2 

E12 = E13 = S23 = 
7 e (4/3, 2) 

Stable for 7 < 4/3 + 6c% + 2c?.. 

3. Equilibrium point 

U=l-9c\- 3cl_ E+ = 3c+ S_ = VSc- 

S12 = S13 = S23 = 

f2^ = 3c^ + c^<l/3 
Stable for 7 > 4/3 + 6c\ + 2c?.. 
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1 1 

Figure 1: Stability zones of the equilibrium points for the case D,\ = 0. Below 
the plane the point 1 is stable, above the plane outside the compact region the 
point 2 is stable, inside the compact region the point 3 is stable. 

4. Equilibrium set (the Kasner Universe) 
U = fi^ = q = 2 

S+ + ^— + ^12 + ^13 + ^23 ~ 1 

The results on stability of the equilibrium points are presented graphically 
in Fig.l. The space (7, c_|-,c_) is divided into several zones in which only one 
stable future attractor exists. 

We can see that the only stable point for 7 < 4/3 is the isotropic one, 
independently on c-|_ and c_. It is possible to strengthen this result and prove 
that all the solutions with > and 7 < 4/3 tend to the isotropic solution 
{flfi = is an invariant subset of our phase space, as it can be seen from Eq. 
(16)). We point out that a physically important boundary case 7 = 4/3 is 
included in the conditions of this statement. The key point is the equation 

^;. = [(6-37)s2 + (4-37)C/]l^^ (19) 

which can be obtained from Eqs. (14) (16). Note, that Eq.(19) follows only 
from the Raychaudhuri and constraint equations and we do not use the evolution 
equations for the shear. This means that the consideration presented below does 
not depend on a particular choice of the angular velocity . 
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First, consider the case 7 < 4/3. Then there is a positive number a such 
that 

> a[Y? + U]%. 

From this inequahty it follows that the quantity Q.^ is strictly increasing if 

^ti{To) > 0. Using this fact and the constraint T,"^ + ilf^ + U ~ 1, one can easily 
understand that the sum E-^ + U, in turn, strictly decreases. As decreasing and 
bounded below by zero, this sum has some limit /3 > 0. Assume that /3 is not 
equal to zero. Then, the quantity il^ goes to infinity: lnf2^(r) — Inrj^(To) > 
a(3{T — To) which contradicts the inequality ^^/^(t) < 1. Thus, for any initial 
conditions satisfying n^(ro) > the solution tends to the point S+ = S_ = 

E12 = 5^13 = ^23 = U = 0. 

In the boundary case 7 = 4/3 it is necessary to be more careful. For 7 = 4/3 
we have 

fl'^ = 2Y?Q.f,. (20) 

So, Q.^ is monotone increasing and, then, tends to some limit value ^um- As- 
sume that fiiim < 1- Using the constraint (14) and the evolution equations for 
U and f2^, we obtain the evolution equation for Y?: 

= 6(ci + C2)UY+ + 2V3(ci - C2)C/S_ + 2S2(S2 - 1). 

OT 

In sequel, we shell use only the fact that the absolute value of the LHS of 
this equation is bounded above: |^S^(r)| < C for all r. Now assume that 
does not tend to zero. Then there is a sequence r„ going to infinity such 
that E^(r„) > c for some c > 0. But the boundedness of the shear square 
derivative means that there is a positive b which does not depend on n such 
that E^(r) > c/2 for t e (t„ — 6,t„ + b) and all n. In this case from equation 
(20) it follows that the quantity goes to infinity: 

r 

Infi^(r) -lnn^(To) = f 2i:'^dt > ^ 2cb. 

TO n:Tn<T 

This fact contradicts the inequality ^^/^(t) < 1. Thus, tends to zero. By 
the constraint 1 = + 0,^ + [/, this means that the quantity U converges to 
1 — riiim > 0. But it is impossible because the point E^ = and U = const > 
is not an equilibrium point of our system for nonzero ci and 02- Thus, 0^ goes 
to 1 and both U and E^ decrease to zero. 

For the brane problem this means that if Vf^u satisfies the ansatz (9) with 
U > 0, then an arbitrary anisotropy of the whole 5-dimensional metric does not 
destroy the future isotropic regime unless 7 > 4/3. 

It is interesting that in the case of 7 = 4/3 without an anisotropic stress 
we have a one-dimensional set of equilibria instead of the isotropic equilibrium 
point. Points of this set have E^ = and are marked by the constant ratio U/^^. 
A nonzero V^i, makes all these points except for U = 0, il^ = 1 unstable. As a 
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result, in the presense of a nonlocal anisotropic stress the "dark energy" U can be 
washed out during the cosmological evolution even in the radiation-dominated 
epoch. 

We finish this section by one remark which is not directly connected with 
the brane cosmology. As we consider the case oi ^\ = 0, the system of equa- 
tions (13) - (16) with Barrow-Maartens ansatz for V^y has the same form (with 
appropriate rescaling of source terms) as the system describing the evolution 
of the Universe in the presence of an anisotropic matter in the ordinary Gen- 
eral Relativity. Since our result on isotropisation does not depend on angular 
velocity Vl" (these values should be determined separately for each particular 
anisotropic source), we can claim that in the presence of an anisotropic matter 
with an anisotropic stress t:^^ proportional to its energy density e and an av- 
erage pressure p equals to e/3 all the models isotropize if 7 < 4/3. Moreover, 
we can consider an anisotropic matter with an average pressure p — {jeff — l)e 
and, using the same procedure, claim that all such models with 7 < -feff tend 
to the isotropic attractor. This fact was previously discovered in the particular 
case of the homogeneous magnetic field [|l^ and for the general source in the 
limit of a small anisotropy |l^ (the "criticality condition" of Ref. |16 is exactly 
7 = Jeff in our notation). 



4 A pure brane regime 

In this section we neglecte the standard matter on a brane. This case describes 
a purely nonstandard brane dynamics when the matter density on the brane is 
large in comparison with the brane tension. The past asymptotic states of these 
models correspond to a brane dynamics near a singularity, future asymptotic 
states describe some intermediate regime of the brane Universe when a nonstan- 
dard asymptotic behavior have already been established but the matter density 
is still larger than the brane tension. 

1'. The isotropic Universe 

U = S+ = E- = S12 = S13 = S23 = 

q = 3j-l nx = l 0<7s^2 

Stable in the future for 7 < 2/3, stable in the past for 7 > 1. 
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c+(37-2) 
Zc\ + 

. 37-2 
94 + 3c^ 

7 e (2/3,1) 

Stable in the future for 7 < 2/3 + 2,c\ + c^_. 
3'. Equilibrium point 

U =l-Qc\- 3c?_ E+ = 3c+ E_ = \/3c_ 

E12 = S13 = E23 = 

f2A=0 3c^+c?_<l/3 
Stable in the future for 7 > 2/3 + 3c^ + c?.. 
4'. Equilibrium set (the Kasner Universe) 
[7 = f2A=0 g = 2 

+ E_ + S;^2 + ^13 + ^23 ~ 1 

Stable in the future for 3c^ + c?. > 1/3. 

Zones of stability of these equilibrium points are plotted in Fig. 2. For 7 < 1 

we can see a complete analogy with Fig.l if 7 is divided by 2. It is natural, 
because the nonstandard matter behaves effectively as a "normal" matter with 
7;vs = 27. Then, to receive expressions for equilibrium points and stability 
zones of the system (13) - (17) with = it is enough to take the correspond- 
ing expressions for fi^ = and replace 7 by 27. In particular, the isotropic 
attractor is stable for 7 < 2/3. 

The part 7 > 1 in Fig. 2 is new (in Fig.l this part would correspond to the 
unphysical region 7 > 2). In this part of the space (7, c+, c_) the stability does 
not depend on 7. As a remarkable feature we can note appearance of the Kasner 
solution as a future attractor for 3c^ + c?_ > 1/3. For 3c^ +(?_ < 1/3 the point 
3' is stable. 

As in the previous section, using the expression 

= [(6-67)S2 + 2(2-37)[/]f}A (21) 



2 . Equilibrium point 

^^0-7X37^ 

94 + 3c2_ ^ 
E12 = Si3 = E23 = 
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Figure 2: Stability zones of equilibrium points for the case = 0. Below the 

solid plane the point 1' is stable, between solid and dashed planes outside the 
compact region the point 2' is stable, inside the compact region the point 3' is 
stable, above the dashed plane outside the compact region the point 4' is stable. 
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we can conclude that all the models with < f^A < 1 and 76 (0, 2/3] isotropize. 

Anisotropic brane attractors may naturally appear in the inflationary brane 
scenario if the transition from nonstandard dynamics to standard one occurs 
at some moment after the inflation. During the nonstandard brane phase the 
condition for inflation is 7 < 1/3 ^ , after the end of the inflation the parameter 
7 grows to 4/3 in typical models. In this case the brane Universe evolves con- 
secutively from one stable regime to another one according to a corresponding 
c+=const, c_=const line in Fig. 2. When the pressure becomes positive (7 > 1), 
the Universe finally reaches one of the two possible attractors depending on the 
parameters c+ , c_ and stays in this regime until the matter density on the brane 
drops below the brane tension and the transition to the standard regime occurs. 

Near the initial singularity the nonstandard matter on the brane dominates 
and we can put the standard matter density flfj. = 0. As we now consider the 
backward time direction, the condition for the stability of the equilibrium points 
becomes A; > 0. The results are follows. 

For the Kasner set two nonzero eigenvalues are Ai = 6 — 67 and A2 = 
1 — 6S+C+ — 2\/3S-C_. Four other eigenvalues are equal to zero. This indicates 
that we have a four-dimensional set of equilibria. The second nonzero eigenvalue 
is positive on some subset of the Kasner set and this subset exists for arbitrary 
(c+, c_). On the other hand, Ai is positive only for 7 < 1. So, the condition for 
stability of the Kasner solution in the past is 7 < 1 independently on (c+ , c_ ) . 

In the opposite case 7 > 1 (a matter with a positive pressure) a stable 
solution in the past direction is the isotropic Universe. This result appears only 
in the brane scenario. In the case of Vf^, — the isotropic character of a brane 
cosmological singularity in the presence of matter with a positive pressure has 
been found earlier by several authors |^ . 

5 Conclusions 

We have investigated the Bianchi I cosmological dynamics on a brane in the 
presence of a nonlocal anisotropic stress Vf_iu satisfying the Barrow-Maartens 
ansatz (9) with U > 0. The dynamics appears to be more complicated in 
comparison with the previously known case of V^i, = 0. However, we can 
mention two important regimes not modified by the anisotropic stress of the 
form studied in this paper: 

• IfO<7<4/3 the anisotropic stress does not prevent the Bianchi I model 
from isotropization. Using the existence of a monotonic function (19) we 
can prove that all such models isotropize. 

• The anisotropic stress does not alter the past asymptotic state which is 
the isotropic Universe for 7 > 1 and the Kasner solution for 7 < 1. 

We also describe the set of future asymptotic states in the case 7 > 4/3 
and possible transient brane dynamics which can exist before the linear energy- 
momentum contribution in the effective Einstein equations becomes important. 
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An important problem for the future development is to consider the U < 
case. 
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